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Abstract 

We derive the low-energy expansion of iZaf' and {ZaY" terms of the polarization operator in the 
Coulomb field. Physical applications such as the low-energy Delbriick scattering and "magnetic 
loop" contribution to the g factor of the bound electron are considered. 
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^ . 1. Introduction 

^, 

One of the predictions of the quantum field theory is a vacuum polarization by an external 
field. An important case thoroughly studied both experimentally and theoretically is the vacuum 
!>■ I polarization effects in atomic field. Methods used for the study of this effect essentially depend 

on the nuclear charge Z\e\. At low Z, the perturbation theory with respect to Za is applicable 
{a - e^ - 1/137 is the fine structure constant, h - c - \). At high Z, the interaction with 
^vq ■ the external field should be taken into account exactly, which can be done with the help of the 

I • I electron Green function in this field. This approach often requires quite involved numerical 

^— ^ i calculations, which usually fail to give the results for low Z. Thus, the two approaches tend to be 

QQ ' complementary. Usually, the perturbative calculations of vacuum polarization eff'ect are limited 

<^^ i by the leading order, since the first nonvanishing correction involves two more loops. Nowadays, 

^ ' the modern methods of calculation of the multiloop integrals are sufficiently powerful for the 

calculation of higher orders in Za. It provides a possibility to compare the results of these 
approaches. 

One of the basic nonlinear QED processes in the atomic field is the Delbriick scattering uJ], 
the scattering of the photon in the Coulomb field due to the vacuum polarization. The amplitude 
of this process in the Born approximation has been obtained long ago for arbitrary energies in 
Ref. ||2|]. At high energies and small scattering angles, when the quasiclassical approximation is 
valid, the amplitude is known exactly in Za, see Refs. 031 |4J, |5l la 13, l8|] . Recently in Ref. {%, 
the Delbriick amplitude has been calculated numerically exactly in the parameter Za at low 
energies. It was shown that the contribution of the higher orders (Coulomb corrections) to the 
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amplitude can be well fitted by the polynomial C4(Za)'^ + C(,(Za)^. The calculation of (Zc)^ term 
in perturbation theory would provide the independent check of the result of Ref. |19[]. 

In present paper, we consider the polarization operator 11'"' {a>, k, q) in the Coulomb field for 
small external momenta 

w ~ |k| ~ Iql ~ Am, (1) 

where A is a dimensionless small parameter, 

A<^1. (2) 

We calculate the expansion of the polarization operator in the Coulomb field in A and Za up to the 
order A'^ (Zdf . The low-energy Delbriick scattering amplitude is readily expressed in terms of 
this operator. This polarization operator is also an essential ingredient of calculations of different 
physical observables in atoms, like Lamb shift and magnetic moment of the bound particle. 
The polarization operator in the external Coulomb field is determined as follows: 

n"" (w, k, q) = \me^ f t/xc^y £/fe-'"'+''"'-''i^ (vac| T/" (f , x) F (0, y) |vac) , (3) 

where J^ - (/ry'^i/r is the electron current and the state |vac) corresponds to the vacuum state in 
the presence of the Coulomb potential Z \e\ jr. In e^ order, we have 

/dp 
—dxdy e'^^-^^Tx [/G (x, y|e - w) y^G (y, x|e)] , (4) 

In 

where G (y, x|e) is the Green function of the electron in the Coulomb field. Due to the gauge 
invariance, the polarization operator obeys the constraints 

k^YV" (w, k, q) = q^IV" (w, k, q) = 0, (5) 

where A;" - ^° - lo. Using these constraints, we can express 11'"' via five independent tensor 
structures, which we choose as follows 



n'"' = /i (g"" k-q- <r) - fT,ef"f^^nae'P''''n 



kp (k - q\ (k - q)^ qj 



■p 



ik - q) 



, kaqp _ ik-q)aik-q)p 

hSafi + /4 2~ ~ J^ 



{k-q? 



(6) 



where n - (1,0) and /] are some scalar functions of w, k and q. 

It is important to note that the low-energy expansion of 11'"' itself and the functions /i_5 is not 
reduced to the multiple Taylor expansion in k' , q' and w, i.e., /i_5 are nonanalytic functions of the 
external momenta. Nevertheless, the expansion in powers of the parameter A from Eq. ^ is still 
possible. Different terms of the expansion come from one of two different regions of integration. 
We separate the contributions of these regions, using the dimensional regularization, in spirit 
of Refs. llldllin . The methods of separation and calculation of the contributions of these two 



regions are given in Sections|2]and[3] respectively. The results and conclusions are presented in 
SectionlH 




Figure 1 : Example of the integral having nonanalytic expansion in q. 



2. Separation of the contributions of soft and hard regions 

In order to demonstrate our method, let us consider the behaviour of the integral 



f 



(X>) 



(q) = J ^^^y(p>A,q) 



J A2 (q + Af (p2 + 1) [(p + q)2 + l] [(p - A)^ + l] ' 

at q^ <s 1. The corresponding diagram is depicted in Fig.[T] The small-q^ expansion of this 
integral, obtained in Ref. [12], has the form: 



/^' (q^) = 2 C„ m (-q^)" + (q^)-^ ^ ^« ^^^ M" 



n=0 



«=0 



(8) 



C„ (D) 



r(n + 3-D/2) 
■ (n + 3 - D/2) (D - 3) 
r(« + 2)r(2-£)/2) 



r(D/2-l) 



r(n + 5-D)r(n + 3-D/2) 
r(n + £)/2)r(2n + 6-D) 



r(2n + 3) 



D„ (D) = r(2 - D/2)r(D/2 - 1> 



2r(« + i)r(« + 3-£)/2) 



(9) 



(10) 



r(£)-2)r(2« + 3) 

The representation dHJ determines 7^®' (q^j at < q^ < 4. The limit q^ — » essentially depends 
on !D. For D > 4, this limit is equal to Co (£)) and can be considered as the value of the function 
J^^^ at q^ - 0. For £) < 4, there is no finite q^ — > limit of the representation ([8]). However, 
the value 7*®^ (0) defined via the analytic continuation with respect to D from the region D > 4 
is still Co (£)) . In other words, the limit q^ ^ is not commuting with the analytic continuation 
with respect to D. The same claim is valid for the derivatives of 7^®' (q^j with respect to q^. 

As a consequence, the integrand expanded in q gives only the terms oc (q^ after the integration 
within the dimensional regularization. These terms correspond to the first sum in the right-hand 
side of Eq. (O: 



/ 



d^Ad^p 



rrD 



«=0 



q''...q'" d"j{p,A,q) 



dq'^ . . . dq'" 



q=0 



I.c„m[-\ 



n=0 



(11) 



Note that the expansion of the massless propagators is valid only in the region A :» q (hard re- 
gion). In order to obtain the rest terms of the expansion (O one has to determine the contribution 
of the region A ~ q (soft region). To separate these contributions, we use the following trick. In 
the soft region the massive propagators can be expanded in both q and A. Let us truncate this 
expansion at some fixed order A^ in q and A: 



;soft(Mp,A,q) = £;*:>,(p,A,q), 



(12) 



7ll(P,A,q) = 



iiWt" 
Now we can identically rewrite J^^^ (q^ j as 

■d^Ad^p 



T=0 



SOth^tC = 0{q"-'). 



rd^Ad^p , , Cd^Ad^p 

= I — -^ — [;(p,A,q)- ;;oft(A^,p,A,q)]+ I — — — ;Wt (A^, P, A, q) . 



(13) 



The contribution of the soft region in the first term in Eq. (fT3] l is suppressed as q^^'^ ^, thus the 
integral of the difference /' (p, A, q) - y'soft {N, p, A, q) is determined by the region A ~ 1 » (7 up 
to O (g'2L(©+A'-3)/2j j (gj-jjj jjj f^j.)-^ (jjg expansion of /soft {N, p, A, q) in q gives scaleless functions 
of A, which vanish after the integration in the dimensional regularization. Finally, we have 



Z 

«=o 



/^£.(P.^.,.^Z/^«.<M.,.. (.4, 



;(«) 



Cd(P'A'q) = 



5«;(p,A,Tq) 



n\dT" 



;l:)(p,A,q)=^"^^(P'^^'^'»> 



n\dT" 



T=0 



Let us now use this approach for the calculation of the low-energy expansion of the polar- 
ization operator in the Coulomb field. The {Zaf^ contribution (A^ is even) is determined by the 
A^-loop diagrams depicted in Fig. |2] In the dimensional regularization, it can be represented in 
the form 



n^;^^j,(a;,k,q) 



/dsi 
(2^ 



dsd^pY\l-,'d''Ai 



'P'^''(fi,p,Ai 



^N-l 



. (D, k, q) . 



Similar to the previous example, there are two different region of integration 

hard region, when s ~ p' ~ A\ ^ j ~ m. 



soft region, when 



s 

A'! 



m. 
Am. 



(15) 

(16) 
(17) 



Again, the expansion of the polarization operator is the sum of the integrals of the expansion 
of the integrand in hard and soft regions. In the coordinate representation, these regions have 
a simple physical meaning. The characteristic size of the electron field fluctuations (the size of 

4 



XA 



A'-l — ^N-2 




Figure 2: The (Za)'^ contribution to tlie polarization operator. Solid lines denote tlie electron propagator, dashed lines 
denote the Coulomb field. 



the electron loop) is of the order 1 jm. Hard region corresponds to the configurations where the 
distance between the Coulomb source and the electron loop is of the order of 1 jm. Soft region 
corresponds to the creation of the virtual electron-positron pair far from the Coulomb source. 
Obviously, there is no contribution from the region where only some of the momenta A„ are hard 
while the rest are soft. In the momentum representation, these regions correspond to massless 
tadpole diagrams which are zero in the dimensional regularization. Thus, the expansion of the 
polarization operator ( fTSl l has the form 

n!': ,„ (w, k, q) = nf; ,„ ., (w, k, q) + If' ^ „ (u, k, q) 



{Zaf 



(Zff)'',hard 



(Zq-)", soft 



^J (2;r)»+' (27r)(^-"^ ^^"^ 4" J (2;r)»+' (27r)(^-"» ^''"' ' ^ 



f,y(n) _ d"^'"' (g, p, Ai , . ■ . Aa,_i , TO), Tk, rq) 



nldr" 



T=0 



qj! 



soft 



„) d"T^-2<^i^y (e, p, tAi , . . . tAa,_i , ™, rk, rq) 



nldr" 



(19) 
(20) 



Note that the simple power counting allows one to estimate the leading terms of the hard and soft 
contribution as 

n^:_.„,. ,(w,k,q)~i^ 



(Za)",hard ' 
(Z»)",soft' 



^(w,k,q)~/^-"*®-'>+'. (21) 

Using Eq. (fTsT l. one can calculate the contributions of the hard and soft regions separately. 



3. Method of calculation 

The contribution of the soft region can be graphically represented as the tree diagram shown 
in Fig.[3] The local multiphoton vertex depicted as a thick dot corresponds to the expansion of 




^ + 

Figure 3: Graphical representation of the soft contribution. 



the fermionic loop with respect to the soft momenta A,, k, q. The expansion is expressed via the 
integrals of the following form 



/ 



ded^p I ^ 2 2 -nr" ii-lYTin-DI2-ll2) 

jtV e -p -wr + ;0 = Ti^rrrr', ■ (22) 

(2;r)»+i ^ ^ ' (47r)(®+"/2p(„)„^2«-B-i 



The remaining integrals over A,- can be easily evaluated in the coordinate representation. Nat- 
urally, the contribution of the soft region can be calculated also with the help of the derivative 
expansion of the one-loop effective QED action. 

The contribution of the hard region determined by Eq. (fTSl ) is expressed in terms of the A^- 
loop tadpoles. In particular, in {Zaf order, the basic integral has the following form 



"J ~(2^ 



dsd^pd^A 1 

(27r)2»+i (e2 _ p2 _ ,„2 + ,-0)". (e2 _ (p_^f _ ^2 + ^-q)"^ (a^)"^ ' 



(23) 



After the Wick rotation s — > ie and rescaling p — » Ve^ + rrfip, A — > Ve^ + ifP-\, we integrate 
over s and obtain 



mi-2yr(y-l/2) Cd^pd^A 
I = i ' 

2 



T(y-l/2) r d^pd^A (-1)"'""^ 

V^r(r) J (27r)2^ (p2 + l)«'((p-A)2 + l)"^(A2)"^' 

y — ni + n2 + riT, - T). 

The remaining integral is the two-loop tadpole in £) = 3 H- e dimensions which can be easily 
expressed in terms of F-functions. 

Performing the similar integration over e in {Zaf order, we express the contribution of the 
hard region in terms of the integrals of the topology (and its subtopologies) depicted in Fig. HI 
The general form of such integral is 

d'°hd'°kid'°hd'°k^D-'">D-""' 



where 



f d^kid^k2d^kj,d^k4D"''D"'" 



Dl,...,A=k\ 4+1, (26) 



^I,...,4 

£•5 = (^1 - k2f , £»6 = {k2 - kif , 
Dj^iki-k^f, Ds = ik4-ki)\ 

Dg = (^1 - k^f , Dio = (k2 - k^f . 




Figure 4: Topology of the integrals required for the calculation of the low-energy expansion of the polarization operator 
in (Za) order. 







(a) icakc (b) ilnflnity (c) iMelon (d) ilumblcr 

Figure 5: Master integrals of the topology, shown in Fig.|4] 




(e) y'cio 



The IBP reduction procedure lll3[ 11411 allows one to express any vacuum integral of the 
considered topology via the five master integrals shown in Fig.|5] On these diagrams the solid 

and dashed lines denote the massive (k^ + l) and massless (k^j propagators. For each loop 
momentum the integration measure is taken as 

Four of these integrals are trivially expressed in terms of F-functions. Their explicit forms are 
presented in Appendix. The only nontrivial master integral is 7cake- After the IBP reduction, the 
master integral 7cake enters the polarization operator with the coeflicient, having the first-order 
pole in the point D - 3. Therefore, we need to determine the OU^j and Ofe') terms of the 
expansion of 7cake- 

We find it convenient to use the recurrence relation with respect to space-time dimension, see 
Ref. II15II . First, we use the Feynman parameterization to obtain the relation 



fiD-2) 'iff j(D-2) 
■'Cake ~ ■'11 nil 1100 



•47' 



1112122200 

<®> +4/ 

1 122121200 ^^■'1212112200 



°-'lll2221200 "^ °-'ll22211200 "^ °-' 1222111200 



t(0) 



AO) 



jiO) 



^■'1212121200 ^ ■'22221 1 1100- 



(27) 



Then, using the IBP identities, we express the integrals in the right-hand side of this relation via 
the five master integrals from Fig.|5] In particular, the last term in Eq. (|27] | can be expressed via 



the master integrals as follows 

,(5+£) dg (5+e) _ ,(5+e) QClovef ,(5+e) 

■'2 - ■'2222111100 ~ "Cake-Zcat-j. "I" ^3 ■'clover 

, r(5+<;) , QTumblei- ,(5+e) QMelon ,(5+e) ^^qn 

+ ainfinity-Zinflnity + ^2 ■'Tumbler ^ ^2 ■'Melon ' ^'^°' 

The coefficients a, are presented in the Appendix. They are chosen to be finite in the limit e — > 0. 
After the reduction, we have the following recurrence relation: 

f(3+f) _ ^y, r(5+f) , ^Clover (5+f) 

■'Cake ~ •"''Cake-'ca^je "I" ^j ■'Clover 

^Tumbler r(5+e) OMelon 



J. h /(5+^> 4- '^Tumbler ,(5+£) '^Melon ,(5+f) .^q^ 

-I- ^Infinity -/infinity "I" ■'xumbler "^ ■'Melon ' '^^^'' 



Again, the coefficients bj are chosen to be finite in the limit e — > and are presented in the 
Appendix. Now we use t 
into Eq. (|27] |. We obtain 



Appendix. Now we use the following trick. Let us express Jq^^^ from Eq. (|28] | and substitute 



_ e^Cake ,(5+e) ^Cake /^Clover ^Clover | ,(: 
aCake e^ \ OCake OCake / 



/ ^Infinity £ ^Infinity \ ,(5+e) ^Cake | ^Tumbler ^Tumbler | ,(5+e) 

^Cake Ibuelon ^Melon | ,(5+e) ,-,f^^ 

"*" ~~r~ ~hZ, n^ , Melon ' '^•'^-' 

e \ OCake OCake / 

Since the integral J^ is finite in D - 5 and the coefficient in front of this integral in Eq. (l30l l 
contains e factor, the first term in the right-hand side of Eq. (|30] | does not contribute in e" order. 
Expanding the coefficient a,, /?, and the four simple master integrals, we obtain 



■^Cake - 6 + ^ 




where C - 0.577... is the Euler constant. Note 


term of Jq^^ "for free". 





^)-.^-j7f].0(.^), (31) 



In order to calculate the O (e) term, let us consider the general solution of the recurrence ( l27l i. 
Taking into account the explicit form of the coefficient /^cake, we obtain 



j?l = Jom 



'Cake 



4 CO (2)+2n) 

tf^./o(© + 2«) ' 



(32) 



where 



, ,^. _ 2^^r(l-£)/2)r(3£)/2-ll/2) 

■^0 (») r(©-2)r(©-3)2 ' ^^^^ 

_ ^Clover _ ^Infinity 

*^Clover ~ ~TJ ! ^Infinity ~ T ? 

f t^Cake ft'Cake 

_ ^Tumbler ^Melon 

*^Tumbler ~ ~~^j ? *^Meloii 



e^ocake e-^o. 



'Cake 



and P (£)) - f (£) + 2) is a periodic function of £). Note that, using the explicit form of the 
simple master integrals, Eq. ( l6Tl l. the sums in Eq. ( |32] | can be checked to converge rapidly. In 
order to fix the function P (£)), we have to calculate the leading asymptotic of ^^,^^1^^. at D — > +00. 
However, the calculation of this asymptotic is not a simple problem. Instead, we may proceed 
in the alternative way by applying the method of difference equations described in Ref. lll6r . 

('f)\ 

According to this method, we derive the recurrence relation in x for -/j-mniioo- 

xiniinoo ~ '"'^■^^-'■'.t+i, 111111100 "*" ■^(■^)' ^''^> 

x(x+ll-3D) 
^^^^ ~ {x + 5- 30/2) (x + 3 - D) ' ^^^^ 

where Fix) is expressed in terms of finite sums of F- functions. The solution of this relation is 



j(fi) MD) ( . 

■'vl 11 11 1100 ^ '-^^ 



F(y) 



y=x -^ 



(36) 



,(2j) rix + 5-3D/2)r{x + 3-D) 



r (x) r (jc + 1 1 - 3D) 

where R'-^^ (x) is a periodic function of x, which can be determined through the J^^ 1 j 1 1 1 100 asymp- 
totic behaviour in x — > 00 limit. At large x, the x-dependence of -/j-imijioo factorizes into 
J d^ki (k^ + lj . Using the large-x behaviour of the function G*®' (x), we find 



t(2)) 

■'.(111 11 1100 



x'-^. (38) 



G(») (x) 

Therefore, for D > 3 we have R(x) - Q and we can use Eq. ( 1361 ) to estimate numerically the 
function P(!D) in Eq. ( l32l i. It should be noted that the recurrence relation ( l34b can be hardly 
applied near 25 = 3 due to the slow convergence of the sum in the right-hand side of Eq. ( |36] |. 
Performing the estimation of P (D) for several non-integer values of D, we find that in all cases 
the value of the function f(D) is compatible with zero up to lO"'". Thus, our ansatz is f (D) - 0, 
and we have 

^^if - ^'~S^^^^ Y.Ti2n.e), (39) 



-r(^)^ 



,3 



T{v) 



8~''-5(v+ I)(l4v3 -h40v2 -h35v+ 10)r(-f - 5^ r(v> 

7r^/2 (6417v5 + 13266v4 + 8368v^ + 980v2 - 585v - 126) r(v) 
6144 v(3v - l)(3v + l)r(-2v)r(3v + 4) cos (f ) sin2(;rv) cos (^) 
71^/2 ^585^3 + 372v2 + 7v _ 12) Y{v)T[^^ tan [^) 

12288(3v - l)(3v 4- l)r(2v + 2) cos (f ) sin {2nv) 
n^l^ (3897v4 + 3ilQv^ + 191v^ - 198v - 54) r(v)2r (^)^ tan [^) 
36864(3v- l)(3v-i- l)r(2vH-2)r(|vH- |)sin(27rv)sin(f ) 
9 



(40) 



The series in this representation are converging rapidly. The first term in T (v), the most slowly 
decreasing one, behaves as 2>^^". Thus, roughly speaking, each two consecutive terms in the sum 
give three more decimal digits of precision. We claim this expression to be valid for arbitrary e. 
For example, one can easily reproduce terms of the expansion of /cake near D - A obtained in 

Using Eqs. (|28|, (l39|, we obtain 



1 *" 

Jf = -—n" (63;r2 - 488 - 72^3) + J] ^2(«) ' (41) 

Tiin) = 



81 

M=l 

(-1)" (56«^ + 80«2 + 35« + 5)n\{n - \)\f 



18(2n+ l)2(3n + l)! 
(_l)"2-4"-i (2340n^ + 744«2 + 7« - b)n\{2n)\)^ 

9n(36n2- l)(«!)2(4n + 1)! 
I6n\2n- 1)!(4«- 1)! 



27(36n2- l)''(6« + 3)! 
X (ll088576n'' + 7641216n^ + 691632n'* - 424512«^ - 79356«^ + 616« + 585) 
2i-V)"7T\i2n)\f{3n)\ 

27n2 (36«2 - \f («!)3(4« + 2)!(6« + 1)! 
X (4489344n'' + 2794176n^ + 206352n4 - 146880«^ - 26032n2 - 28n + 153) 

AnH2n)\{An)\ (//2„-i + 2//4„-i - 3//6„+3) 
27n2(36n2- l)(6n + 3)! 

X (l02672n^ + 106128n'* + 33472«^ + \9(i0n^ - 585n - 63) 

(-l)";r^((2«)!)^(3«)! (3//„-i - 10//2„-i - 3//3„+i +4//4„+2 + 6//6«+3) 
54n2(36n2_ 1)(„!)3(4„ + i)!(6n + 1)! 

X (31176«'' + 15480n^ + 1594«2 - 198« - 27) , 



'2 
-1 



£^ = 0.0516516357945 .... (42) 



Here H„ - ^11=1 ^ i^ ^ harmonic number. Unfortunately, we have not been able to express the 
sums in Eq. (HTt in terms of ^-functions and alike. However, the numerical convergence of the 
above series is perfect and we obtain from Eqs. (ISTT l. (HTI) 



^4 



4lke' " ^ ~ ^ ^ 58.3184377060 . . . + O (e^) . (43) 

Methods of calculation and values of multiloop vacuum integrals in arbitrary space-time di- 
mension D are of independent interest, because these integrals appear as parts of the amplitudes 
for various physical processes: from QCD and QED radiative corrections lilS1 to the thermody- 
namics of finite temperature QCD-like theories il9tl . In particular, the master integrals in Fig.|5] 
enter the basis intensively used in modern QCD calculations Il20ll2lll22ll23ll . 
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4. Results and Conclusion 

The perturbative expansion of the form factors in Eq. (|6]l has the form 

fi^ J] a(Zarf^"\ 



(44) 



n=2,4,- 



We have calculated the low-energy expansion of /j and /j up to (9 (/l^j, which corresponds to 
the expansion of the polarization operator up to O (/!'* j. Using Eq. ( fTsT l. we obtain in (Za)^ order 



,(2) ^ 



7 ^ 1 |k-q| ^ 



1 



f(2) 



2(4!)2 6! m 2H5l)^m^ 

73 ll|k-q| 



liZ^2+49k.q-^(k2+q^) 



22(4!)2 

1 

25(5!)2m2 



6! m 
10923 



w^ + 3095k- q+ 1111 (k2 + q2) 



/, 



{2) 



7 |k -q| 587 (k - q)^ 



/ 



6! m ■ 23(5 !)2 m^ 
573 a/ 



<2) 

4 24(5!)2m2 

<2) 



r^^; 

•^5 AI 



4 |k-q| 1369 (k-q)^ 



m 



25(5 !)2 



(45) 

(46) 
(47) 
(48) 
(49) 



The O(A^) terms of f^^ have been obtained in Ref.(|2|]. When k = w = 0, the expansion of 

/j (q) agrees with the exact result obtained in Ref. i24ll . The next-to-leading terms, proportional 
to |k - ql, come from the soft region and exhibit the above-mentioned nonanalytic behavior. 
In (Za)'* order we obtain 



,(4)_1 /2,(5) 4 5 13 

^' -2^b-^2 -3^3 + ^6 + Y^ 



(50) 



1 



109 

'2^ 



yf + 



475 



.^3. 



4111 



■^2 + 



2575 \ 



72 
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Figure 6: The "magnetic loop" contribution to tlie bound electron g factor. 
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In {Za)'^ order, the nonanalytic contribution to the form factors is suppressed as 0{A^) and thus 
is far beyond the accuracy chosen. Note, that the technique used in this paper can be applied 
without modification to the calculation of the higher terms of the low-energy expansion in {Zaf' 
and {Za)'^ orders. 

The Coulomb corrections to the form factors f\ and /2 in order /l" were calculated in Ref. (|9|] 
numerically. Although the interaction with the Coulomb field was taken into account exactly, it 
turned out that the results can be well fitted by the polynomial function of Za: 

/i = 5(^(Z«)' + 3-35 ■ \Q-\Zaf + 1.6 ■ \Q-\Zaf, (55) 

fi = -^^{Zaf - 3.55 ■ \Q-\Zaf - 2.1 ■ \Q-\Zaf . (56) 

In O (/I'M order, our results ( fSOl l. ( BTI ) for the (Za')'*-order corrections numerically coincide with 
those of Eqs. ( fSSl ), ( |56] | with an accuracy of a few percent. 

As the demonstration of possible applications of our result, let us calculate the O ({Zaf\ 
contribution of the "magnetic loop" to the g factor of the bound electron, see Fig.|6l The corre- 
sponding correction to the g factor of the electron in nLj state has the form 12511 

— = -7^ 7T — /i iq) G (q) , 

go n(2K-l)J q 

f\ iq) = /i lk=o;=o ' G{q)^ I dra (r) b (r) (sin qr - qr cos qr) , 

go = r-^, ^ = (7+ 1/2) sign (L- 7), 
2k + I 
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where a and b are determined by the form of the bound electron wave function 

"^^ ' \ ib(r)njLM(n) 

The characteristic scale of the function G (q) is mZa/n <s: m, so we have two regions of integra- 
tion: q ~ mZa/n and q ~ m. In the leading order, only the first region is essential, the contribu- 
tion of this region is of the order O UZa)^). The leading correction ~ O {(Za)^] for L ^ states 
also comes from the region q ~ mZcc/n, while for L = the whole interval mZa/n ^ q ^ m 
is essential. This correction has been found in Ref. ll25ll . Note that the integrals in Eq. (18) of 
Ref. |25] can be taken analytically, and the correction to g factor up to the order O ((Za)^) can be 
represented as 



80 /(Zaf+lZaf 



la iZaf 



288nV(y+l)(27+l) 



+ 5l=q 



+ Su=o 



AaiZaf 1 
1357r n^ 



log 



2Za 



la{Zaf 



641 («H-1)(4«-1) 

240 " ""^ 6^?2 



1 



\5Tin\lL+X)(lK- 1)2 \L(Lh- 1) 



(57) 



In order to find the next-to-leading correction, we separate the contributions of the two regions 
similar to what has been described above. The details of this calculation will be presented else- 
where. It turns out that the complete 0\{Za)'\ result for the correction to g factor can be ex- 
pressed via several first term of expansion of the function f\ near q-^, namely 
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Now, owing to Eq. (ISOl l, we have the last essential ingredient to obtain the correction. Using 
Eqs. (gill, dlOll, we obtain 
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In particular, for IS 1/2 and 2Pi/2 states we have 



— = 1.62xlO"^Q'(ZQfr +9.4314xlO""Q'(ZQ'r In 2.67 

go I IS \ 2Zq; 

+ 4.1 X IQ-^aiZaf , 

— I = 1.62 X 10"^q'(Zq')^ +5.8946 X lO^^aCZa)* 
golip 

+ 3.26xlO"^a(Za)^ (60) 

The contribution of the oUZay) term is rather essential, e.g., for Z - 6 (carbon) the ratio of 

this term to O ((Za)^) term for 15 1/2 state is -0.81. The last term in Eq. ( |59] l corresponds to the 
contribution of the electron loop with four Coulomb exchanges. It is interesting to compare the 
magnitude of this term with that of the first two terms. As it was claimed in Ref. ll25ll this term 
appears to be numerically small. E.g., for the ground state, the contribution of the last term is 
only 2.2 percent. 

A. Appendix 

The explicit form of the four simple master integrals from Fig.|5]reads: 
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The coefficients in Eq. ( I28i l are 
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(62) 



Here x„ - F (x + n) /F (x) = x(x + I) . . .(x + n - 1). The coefficients in Eq. (|27] | are 
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